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ABSTRACT 



Context. As a result of the numerous missions dedicated to the detection of Gamma-ray bursts (GRBs), the observed properties of 
these events are now well known. However, studying their parameters in the source frame is not simple since it requires having 
measurements of both the bursts' parameters and of their distances. 

Aims. Taking advantage of the forthcoming Catalog of the High Energy Transient Explorer 2 (HETE-2) mission, the aim of this 
paper is to evaluate the main properties of HETE-2 GRBs - the peak energy [Ep^ai^), the duration {Tgo) and the isotropic energy (£,,„) 
- in their source frames and to derive their unbiased distribution. 

Methods. We first construct a complete sample containing all the bursts localized by the Wide-Field X-ray Monitor (WXM) on-board 
HETE-2, which are selected with a uniform criterion and whose observed parameters can be constrained. We then derive the intrinsic 
Epeak, Tgo and Eiso distributions using their redshift when it is available, or their pseudo-redshift otherwise. We finally compute the 
'volume of detectability' V,„„v of each GRB, i.e. the volume of the universe in which the burst is bright enough to be part of our 
sample, and the corresponding number of GRB within their visibility volume Nvmax, in order to derive a weight for each detected 
burst accounting both for the detection significance and the star formation history of the universe. 

Results. We obtain unbiased distributions of three intrinsic properties of HETE-2 GRBs: ii™^^., Tg""' and the isotropic energy of 
the burst. These distributions clearly show the predominence of X-ray flashes (XRFs) in the global GRB population. We also 
derive the rate of local GRBs: R^^ > II Gpc^^^yr"', which is intermediate between the local rate obtained by considering only the 
'high-luminosity' bursts (~I Gpc^^yr"') and that obtained by including the 'low-luminosity' bursts (> 200 Gpc^'yr"'). 
Conclusions. This study shows that the XRFs are predominent in the GRB population and are closely linked to the 'classical' GRBs. 
We show that HETE-2 detected no low-luminosity GRB like GRB 980425 or XRF 060218, due to the small size of its detectors, 
excluding this type of burst from our statistical analysis. The comparison of the GRB rate derived in this study with the known rate 
of Type Ib/c supernovae clearly shows that the progenitors of SNe Ib/c must have some special characteristics in order to produce a 
gamma-ray burst or an X-ray flash. 
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1. Introduction 



of X-ray flashes. 



In recent years we have learned that long Gamma-ray bursts 
{Tgo > 2 sec.) are associated with the death of massive stars 
(M > 20-30 Mq). This origin has been clearly established by 
the association of a few nearby GRBs with Type Ib/c supernovae 
(SNe Ib/c): ORB 980425 detected by B eppo-S AX (Pian et al. 
[T999] l and SN 1998bw (Galama et al. IT9981 Kulkarni et al. 
[T998] l, ORB 021211 detected by HETE-2 (Crew et al. l2002l 
|2003T l and SN 20021t (DeUa Valle et al. 2003J, GRB 03 0329 
also detected by HETE-2 (Vanders pek et al. I2U031 l2004b and 
SN 2003dh (e.g. Matheson et a l. l2003l Stanek et al. l2003l 
Hjorth et al. 2003a; Mazzal i et al.|2003l l. GRB 031203 detected 
by INTEGRAL (Gotz et al. i2003l l and SN 20031w (Malesani et 
al. [2004 ). GRB 050525A detected by Swift-BAT (Band et al. 
|2005] l and SN 2005nc (DeUa Valle et al. 2006b). The association 
of flie nearb y GRB 060218 detected by Swift -BAT (Cusumano 
et al. 120061 ) widi SN 2006aj (Masetti et al. l2006l Modjaz et 
al. IJOOe"; Pian et al. 2006; Campana et al. '2006'; Sollerman et 
al. [2006; Mirabal et al. 2006b. Cobb et al. ,2006) provided a 
remarkable example of this paradigm. Even though a massive 
star's explosion can explain the production of both the GRB and 
the supernova (see for instance the coUapsar model proposed 
by Woosley |1993 ; see also Woosley & Bloom 2006; Delia 
Valle |2006al l. these two phenomena are profoundly different in 
nature: the GRB is due to an ultra-relativistic outflow generated 
by a newborn black hole or magnetar, while the supernova is 
powered by the radioactive decay of ''^Ni within a massive shell 
of matter ejected at sub-relativistic velocities (e.g. Soderberg 
et al. I2006d] l. The question of the link between GRBs and SNe 
has recently become more complicated by the discovery of two 
low-redshift, long gamma-ray bursts - GRB 060505 at z =0.089 
and GRB 060614 at z =0.125 - which were not associated with 
a supernova (Gehrels et al. 2006; Fynbo et al. |2006bl Delia 
VaUe et al l2006dl Gal- Yam et al. 2006, ). It is thus fair to say that 
despite very significant advances in our understanding of GRBs, 
the global picture continues to escape us and key issues remain 
to be elucidated. Concerning the connection between GRBs and 
supernovae, the observational clues are based on two facts: the 
clear association of some GRBs with Type Ib/c supernovae, and 
the evidence, based on statistical studies, that the majority of the 
9000 Type Ib/c supernovae exploding each year in the universe 
do not produce GRBs. This last point was addressed by Berger 
et al. ( 120031), Soderber g et a l. (2006a 2006b] I2006cl I2006dl 
|2004bD and Delia Valle ( |2005ll2006a. ,2006c J . 

While the SNe Ib/c population is rather well sampled in 
our local environment, the situation is quite different for GRBs, 
which need a complex series of detections to be fully char- 
acterized: detection of the high-energy signal in space, quick 
transmission to the ground, identification of the afterglow at 
X-ray, optical or radio wavelengths, and measurement of the 
redshift. This situation makes our understanding of the GRB 
population much less secure than for supernovae. Indeed, the 
observations by Beppo-SAX, HETE-2, and now Swift have 
demonstrated the diversity of the GRB population, which is 
composed of classical GRBs, short-duration GRBs, X-ray 
flashes (XRFs), and low-luminosity GRBs. In this classification, 
the 'classical GRBs' appear to be the 'high-luminosity' part 
(HL-GRBs) of a population which could consist mostly of low- 
luminosity or subluminous bursts (LL-GRBs), like GRB 980425 
(Tinney at al. (T998] l and GRB 060218 (Mirabal et al l2006al l. or 
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In this paper we take advantage of the broad energy cov- 
erage of HETE-2 instruments (2-400 keV) to discuss the 
global properties of the GRB population, and to infer a rate 
of GRB+XRF events which can be compared with theoretical 
predictions and with the known rate of Type Ib/c supernovae. 
This work is partly based on the forthcoming catalog of HETE-2 
GRBs (Vanderspek et al. i2008b . For a description of die H ETE-2 
mission and its instrumentation, see Ricker et al. ( I2001I I, Atteia 
et al. (I2003al l and Villasenor et al. (l2003l l. 
Our work relies on the construction of a complete GRB sample 
containing all the long-duration bursts localized by the Wide 
Field X-Ray Monitor (WXM) and having a signal-to-noise ratio 
(SNR) larger than a given threshold in FREGATE (FREnch 
GAmma-ray TElescope) or in the WXM. For each burst in 
this sample, we compute the intrinsic properties at the source 
(Tgo, Epeak, and Eiso) by correcting the observed properties 
for the effect of the redshift. When a spectroscopic redshift is 
not available (62 GRBs out of 82), we use the pseudo-redshift 
following the method of Pelangeon et al., i.e. the estimate of the 
redshift which is partly based on both the Epeak-Eiso (e-g- Amati 
et al. 2002) and Epeak~Liso (Yonctoku et al. ,2004) correlations, 
and derived from the spectral properties of the prompt emission 
(for more details see Pelangeon et al. |2006al and Pelangeon 
& Atteia 2006b). We emphasize that this procedure has little 
impact on the final results since it introduces an additional 
uncertainty which is much smaller than the intrinsic dispersion 
of the parameters under study. In a second step we derive the 
'visibility distance' for each burst (z„,iia), the distance at which 
the SNR of the burst reaches the threshold of our analysis, and 
the number of GRB within its visibility volume (Nvmax)- We 
then attribute to each GRB a weight W = 1 INvmax, accounting 
both for the detection significance and the star formation 
history of the universe. Consequently, this method enables us 
to renormalize the global distribution, taking into account the 
true rate of occurrence of each type of GRB. We finally derive 
the GRB rate detected by HETE-2, and we discuss the relative 
importance of XRFs and classical GRBs in the overall GRB 
population. 

This paper is organized as follows: Section |2] is devoted 
to the description of the HETE-2 GRB sample. Section [3] 
references and describes the tools necessary for our study. 
The three parameters studied here, Ep^ak, Tgo, and E/so, are 
discussed in Sections |4] |5] and |6] respectively. The remaining 
three sections are dedicated to the interpretation of the results. 
In Section|2l we check the existence of correlations between the 
intrinsic parameters. In Section[8] we show that our study gives 
some clues to the nature of X-ray flashes in terms of intrinsic 
energetics and distance-scales. Finally, in Section |9] we discuss 
the rate of GRBs in the local universe and compare it with 
previous estimates obtained by other authors. A summary of the 
main results obtained in this paper is presented in SectionfTOl 

2. Obtaining of a complete sample of HETE-2 bursts 

2.1. Instrumentation 

Since its launch in October 2000, the HETE-2 satellite (Ricker 
et al. 120011 ), dedicated to the detection and observation of 
Gamma-Ray Bursts, has detected 250 events classified as GRBs 
(Vanderspek et al. 2008). Thanks to the combination of its 3 
instruments, the FREnch GAmma-ray TElescope (FREGATE, 
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Atteia et al. |2003al), the Wide-field X-ray Camera (WXM, 
Shirasaki et al. '2003) and the Soft X-ray Camera (SXC, 
Villasenor et al. 2003 ), HETE-2 has a broad energy range cover- 
ing the hard X-rays and y-rays (2-400 keV), allowing the detec- 
tion of both classical and soft GRBs (XRFs). Moreover, HETE- 
2 was designed to provide GRB positions to the community 
through the GRB Coordinates Network (GCN) (Barthelmy et 
al. [2000 ) with an accuracy of one to several arcminutes, within 
seconds of the trigger Consequently, multi-wavelength ground 
follow-up was done for most HETE-2 GRBs, leading to the de- 
termination of 25 spectroscopic redshifts. 

2.2. Burst selection 

It is of prime importance for our study to construct a GRB sam- 
ple for which the detection criteria are fully understood and can 
be reproduced. This is essential for the determination of the visi- 
bility volume of each GRB in the sample. In addition it is neces- 
sary to be able to compute for each burst the three parameters of 
our study: Epeak, 7^9o, and the fluence. This section explains the 
construction of the sampleQ. 

2.2.1 . Availability of the spectral data 

We considered all the 250 events that constitute the total sample 
detected by HETE-2, i.e. both triggered and untriggered bursts, 
thus reducing the bias due to trigger algorithms. Then, we se- 
lected the bursts for which the angle of incidence was measured 
by any experiment. Its knowledge is necessary to perform the 
spectral analysis of the bursts and to derive their £^^0^ and Ejso. 
Of the 250 GRBs, 132 have no accurate localization and were 
rejected, leaving 118 GRBs. 

Among these bursts, we have also rejected GRB 0408 10 because 
this long and intense burst occurred as the instruments were be- 
ing shut down. Consequently only the precursor was detected by 
HETE-2. 

2.2.2. Technical problems 

Of the remaining 117 GRBs, 2 bursts were rejected due to tech- 
nical problems: GRB 010110 was detected and localized by the 
WXM when FREGATE was not working. It is thus impossible 
to constrain the spectral parameters with only the WXM data. 
The second burst rejected is GRB 021 1 13. At the moment of its 
detection there were onboard software problems that did not per- 
mit the trigger to reach the ground. The spectral data are thus not 
available. 

2.2.3. WXM localization 

Considering the 115 remaining GRBs with a known angle, we 
removed 24 GRBs for which the angle was obtained through the 
localization by the Inter Planetary Network (IPN) or the Swift 
satelUte (Gehrels et al. |2004| l. Indeed, as the IPN and Swift have 
different sensitivities, we would have added a bias to our sample 
by considering them. 

Moreover, in order to have bursts well within the field of view 
of the WXM, we have considered an 'incident angle limit' of 
45°. This angle corresponds to a limit of detection for the WXM 
(Shirasaki et al. |2003l l, and to a decrease of the effective area for 



' The order of the steps described in this Section is not important; we 
could have performed them in any order and obtained the same sample. 



FREGATE (Atteia et al. l2003al l. With this cut, only one burst is 
rejected: GRB 020201, whose angle is 55°. This step leaves 90 
GRBs in the sample (see Fig.[T]i- 

2.2.4. Short-duration bursts 

We also rejected the two bursts classified without ambiguity 
as short-duration bursts: GRB 02053 1 (Lamb et al. 2003) and 
GRB 050709 (Villasenor et al. l2005l) . This is based on the fact 
that this class of burst is probably not associated with the same 
progenitors as the long GRBs. As one goal of this paper is to 
discuss the rate of long GRBs in the local universe and to com- 
pare it with the rate of SNe, we have excluded the short-duration 
bursts, leaving 88 bursts in the sample. 

2.2.5. Threshold cut 

Finally, in order to have homogeneous detection criteria for all 
the events of the sample, we removed the GRBs with a signal- 
to-noise ratio (SNR) lower than a given threshold in FREGATE 
and in the WXM. These SNR were computed as follows: 

- we measured the peak count rate of each burst. For 
FREGATE we used two different time resolutions (1.3 and 
5.2 s) and two different energy ranges (6-80 and 30-400 keV) 
corresponding to the bands used by the instrument to trigger. 
For the WXM we used three different time resolutions (1.2, 
4.9 and 9.8 s) and the total energy range of the WXM (2- 
25 keV). 

- we computed a SNR for each burst in all the combinations 
of time resolution and energy range. We kept only the GRBs 
for which the SNR exceeded 7.9 in FREGATE or 5.8 in the 
WXM (in at least one combination of the time and energy 
ranges, see Fig.|2]i. 

This last criterion removed 6 bursts with a low SNR: 
GRB 011103, XRF 020903, GRB 030323, GRB 030706, 
GRB 040131 and GRB 040228. 

Our final sample contains 82 GRBs whose properties are 
given in Tables|2]and[3] 

2.3. The cases of XRF 020903 and GRB 030323 

These two bursts have a spectroscopic redshift but they are re- 
jected from our sample due to their SNR just below the thresh- 
olds of FREGATE and the WXM (Fig. ^. As they have spec- 
ti-o scopic redshifts - z = 0.25 for XRF 020903 (Soderberg et 
al. |2004al l and z = 3.372 for GRB 030323 (Vreeswijk et al. 
I2003I I - we can argue that even if they had been present in our 
sample, the results presented in this paper would actually have 
been reinforced. Indeed, situated at a moderately high redshift, 
GRB 030323 has a large visibility volume and thus a low weight. 
On the other hand, the nearby XRF 020903 has a small visibility 
volume and hence a large weight. Since this is an XRF with a 
low Epeak, adding it to our GRB sample would have increased 
the predominance of X-ray flashes which is shown in Section|4] 
and strengthened the conclusions of this paper 

3. Deriving the intrinsic properties of [-IETE-2 GRBs 

3.1. Burst parameters 

We have used the spectral parameters and durations available 
in the HETE-2 Catalog (Vanderspek et al. l2008i l which summa- 
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Fig. 1. Boresight angle cumulative distribution of the 1 15 GRBs 
localized by HETE-2/WXM, the Interplanetary Network or 
Swift. The dark dashed vertical line represents the 'angle limit' 
of 45°, corresponding to the limit of detection by the WXM. The 
GRBs localized by the WXM are represented with black open 
circles, the ones localized by the IPN are in blue filled triangles, 
and the ones localized by Swift are in dark-red filled stars. 



rizes the information on GRBs detected by HETE-2 during the 
entire mission and thus complements the previous spectral anal- 
yses based on smaller samples of HETE-2 GRBs performed by 
Barraud et al. ( |2003l I2004D and Sakamoto et al. ( |2U05] l. 

3.2. Burst distance 

At the end of the 90s, the growing sample of GRBs with 
measured redshift allowed new types of studies: searches for 
correlations between physical quantities - characterizing the 
gamma-ray bursts' light-curves and/or spectra - in their rest 
frame. Following the discovery of correlations between Eiso, the 
isotropic equivalent energy emitted by the burst, and various in- 
trinsic GRB properties, it has been suggested that some GRB 
observables could be used as luminosity indicators, and hence 
as redshift indicators. The most used are the lag-luminosity 
(riag-Lpjso) correlation (Norris, Marani & Bonnell 2000), the 
variability-luminosity (V-Lpjso) correlation (e.g. Fenimore & 
Ramirez-Ruiz |2000l Reichart & al. 1200 II Lloyd-Ronning, Fryer 
& Ramirez-Ruiz 12002a ) and the Epgak-Ejso correlation - sug- 
gested by several authors (see e.g. Lloyd, Petrosian & Mallozzi 
120001 Atteia 2000) Lloy d-Ronning & Ramirez-Ruiz [2002bD - 
before Amati et al. (120021 1 firmly established it with a study of 12 
GRBs with known redshift. Other correlations were found later, 
such as the Epeak-Lpjso correlation (Yonetoku et al. '2004'), the 
Epeak-Ey correlation (Ghirlanda, Ghisellini & Lazzati 2004 ) and 
the E peak-Lpjso-ToAs correlation (Firmani et al. 2006). Recently, 
in his construction of the Hubble Diagram from a sample of 69 
GRBs with redshifts, Schaefer (120071 1 showed that the use of 
five of these correlations simultaneously - lag-luminosity (r/ag- 
Lpjso, Norris et al. 2000); V -Lpjso (Fenimore & Ramirez-Ruiz 
I2OOOI 1: Epeak-Lpjso (Schaefer l20031 l: Epeak-Ey (Ghirlanda et al. 
120041 1 and the minimum light-curves rise time versus luminosity 



Fig. 2. Signal-to-noise ratio (SNR) cumulative distribution of 
the 88 GRBs localized by the WXM and within 45° of the bore- 
sight. The threshold of FREGATE (SNR=7.9), is represented 
with the right-most vertical magenta dashed line, whereas the 
threshold of the WXM (SNR=5.8) is represented with the left- 
most vertical magenta dashed line. The different symbols show 
the configuration of time resolution and energy range leading to 
the highest SNR for each burst: a given shape is associated with a 
given energy range and a given color corresponds to a given time 
resolution: [WXM/2-25 keV/1.2 s]= light-green filled circles; 
[WXM/2-25 keV/4.9 s]= Ught-green open circles; [WXM/2- 
25 keV/9.8 s]= dark filled circles; [FREG/6-80 keV/1.3 s] = 
red filled triangles; [FREG/6-80 keV/5.2 s]= red open trian- 
gles; [FREG/30-400 keV/1.3 s]= blue filled squai-es; [FREG/30- 
400 keV/5.2 s]= blue open squares. Note that there is only one 
light-green filled circle and one light-green open circle (left part 
of the figure) due to the fact that only one burst has its highest 
SNR in the configuration [WXM/2-25 keV/1.2 s] and one burst 
in the configuration [WXM/2-25 keV/4.9 s]. 



{trt - Lpjso, Schaefer I2OO21Q - led to more reliable luminosity 
distances than the use of only one distance indicator 
The design of HETE-2 permitted the study of the energetics at 
work in GRBs, both in the observer and in the source frames. 
For example, the Ep^ak-Eiso correlation (Amati et al. |20021 l was 
confirmed and extended at lower energy, thanks to the sample 
of HE TE-2 X-ray rich GRBs (XRRs) and X-ray flashes (Amati 
[20031 Lamb, Donaghy & Graziani 2004a; Lamb et al. 20041 
Amati I2OO6I 1, particularly XRF 020903 (Sakamoto et al. 20041. 
The HETE-2 sample is thus appropriate for the determination 
of pseudo-redshifts based on the Epeak-Eiso correlation. A first 
attempt in this direction was proposed by Atteia ( |2003bl l, which 
was later revised by Pelangeon et al. (2006a) to take into account 
the complex lightcurves of some long GRBfl 
The accuracy of this redshift indicator and its use in this study is 
discussed in the following paragraph. 



- See e.g. Ghirlanda, Ghisellini & Firmani i2006l l or Schaefer {2007) 
for a complete summary of these correlations. 
^ See also ht tp: //www, ast . obs-mip. fr/grb/pzj 
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3.3. Justification of tine pseudo-redsiiifts used in this study 



10r 



pscudo-z / z 



As shown in Table |3] 18 GRBs contained in our sample have 
a spectroscopic redshift and 2 have a photometric redshift 
(GRB 020127 and ORB 030115). We computed a pseudo- ^ ^ ★ 
redshift for the remaining 62 (among the 20 bursts having a ^ i 

known redshift, 14 were used for calibrating the pseudo-z). 
In order to test the impact of the use of pseudo-redshifts in this 
work: 




0.1 



- we computed the pseudo-redshifts of the 20 GRBs having a 
spectroscopic or photometric redshift, 

- for these same 20 GRBs we computed the ratio between the 
pseudo-redshift and the redshift, and the ratio of the pseudo- 
intrinsic properties derived with pseudo-redshifts to the in- 
trinsic properties derived using the redshift (Fig.|3]l. 

We note that the dispersion of the ratio between the pseudo- 
redshifts and the redshifts is smaller than a factor of 2, 
except for two outliers: GRB 020819 and GRB 051022 for 
which the spectroscopic redshift differs from the pseudo- 
redshift by factors of 2.95 and 2.15 respectively. For these 
20 GRBs, the dispersion between the luminosity distance 
estimated with the pseudo-redshifts {DiesHm) and the lu- 
minosity distance measured with the redshifts {Di^meas) is: 
CTfl^ = log(DL^e„i„/DL,meas) = 0.125 dcx, i.e a factor of l.Sfl 
The dispersion is also smaller than a factor of 2 for the Epeaii 
(Fig.|3] top-right panel) and the Tgo parameters (Fig.[3j bottom- 
left panel). This is not surprising since these two quantities 
vary as (1 -H z), reducing the impact of the redshift errors. Etso, 
which is more redshift-dependent than Epeak and the duration, 
is more scattered and its dispersion is larger than a factor of 
2. This corresponds to the bursts having the highest difference 
between their spectroscopic redshift and their pseudo-redshift: 
GRB 020819, GRB 0409 12B and GRB 051022 have a pseudo-z 
estimate two to three times higher than their redshift; and 
GRB 020124, GRB 030115 and GRB 050408 for which the 
pseudo-z is lower than the spectroscopic redshift by a factor of 
2 (Fig. [3] bottom-right panel). 

We consider that the good agreement of pseudo-z with spec- 
troscopic redshifts for the 20 GRBs having a measured z is 
encouraging. One may question however if the pseudo-z values 
computed without knowing the spectroscopic redshift are as 
accurate as those computed when the redshift is known - i.e. do 
pseudo-z have predictive power? - Unfortunately, most of the 
bursts for which we have computed a pseudo-z do not have red- 
shift measurements, so comparison with the redshift-estimates 
cannot be done. However, in some cases the pseudo-z values 
we issued in GCN Circulars were followed by a spectroscopic 
measurement. For instance, we determined the pseudo-redshifts 
of GRB 05052^1 and GRB 0701250 before their redshifts 



If we do not consider the two outliers GRB 020819 and 
GRB 051022 the dispersion is croi^ = 0.089 dex, i.e a factor of 1.23. 

' The first GCN Circular containing the computation of a pseudo- 
redshift in nearly real time was done for this burst, using the spectral pa- 
rameters obtained by Konus-WIND (Golenetskii et all20Q5all. We found 
z = 0.36 ± 0.10 (Atteia & Pelangeon 2005b I whereas the spectroscopic 
redshift obtained by spectroscopy of the host galaxy and published the 
day after was z = 0.606 (Foley et al. I2005t . The difference between 
the pseudo-redshift and the redshift was due to an incorrect published 
fluence (Golenetskii et al. I2005bb . As soon as the correct value was 
available, we recomputed the pseudo-redshift of this burst and found 
z = 0.64 ± 0.10 (Atteia & Pelangeon 2005c I. 

* We have computed a redshift-estimate of this bright burst using the 
spectral parameters derived from two different instruments, Rhessi (z = 



pseudo-T""' / T" 

*^ 90B ' 9( 



1— 



0.1^ 



1 

redshift 



1 



1 



1 

redshift 



Fig. 3. Ratio between the pseudo-redshift and the redshift for 
the 20 GRBs contained in our sample for which both quanti- 
ties are available (upper-left panel). The ratio of the 3 pseudo- 
intrinsic and intrinsic parameters studied are also shown, respec- 
tively Epeak (upper-right panel), Tgo (lower-left panel) and Eis„ 
(lower-right panel). In all plots, the filled region corresponds to 
a deviation from the equality between the two parameters by a 
factor of 2, and the error bars are at the Icr level (stat.H-syst.) 



were known. The spectroscopic measurements have confirmed 
our estimates. Moreover, two other GRBs detected in 2003 
and 2004 by HETE-2 had a pseudo-redshift (Pelangeon et al. 
I2006al l but no redshift. Some later observations of the host 
galaxies resulted in the determination of their spectroscopic 
redshifts, in agreement with their pseudo-redshifts: Rau et al. 
(I2005l l measured the redshift of the host galaxy of GRB 03053 1 
to be z = 0.782 + 0.001 (|_= 0.64 + 0.15, Pelangeon et al. 
I2006al l. and Stratta et al. (I2007I I obtained the redshift of the 
X-ray flash 040912: z = 1.563 + 0.001 from the [Oil] line of 
the host-galaxy spectrum, also consistent with its previously 
determined pseudo-redshift (z - 2.90 + 1.60, Pelangeon et al. 
E006al l. 

In order to quantify the effect of the uncertainty from the 
pseudo-redshifts on the intrinsic distributions, we performed 
the following test: we artificially increased the errors on the 
pseudo-z values which are taken into account in the construction 
of the intrinsic parameters distributions (see Section l43T l. This 
extra dispersion takes on three values corresponding to factors 
of 1.4, 2.0 and 3.0. We tested the impact of this extra dispersion 
on the intrinsic Ep^ak distributions ('simple' and 'unbiased'), 
and thus we show that even if the pseudo-z were only accurate 
to a factor of 2 or 3, this would not significantly change the 
intrinsic parameter distributions and all the subsequent results 
obtained in this paper (see Sect.HJl. 



It was recently proposed by Butler et al. ilQQli that the 
correlations between the intrinsic GRB parameters found in 



1.63 ± 0.80, Pelangeon & AH:eia |2007a) and Konus-WIND (2 = 1.34 ± 
0.30, Pelangeon & Atteia 2007b'). Both values were consistent with a 
redshift of ~1.5. The spectroscopic redshift z ~ 1.54 was measured 
several days after our estimate (Fox et al. 120071 see also Cenko et al. 
12008) . 
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the pre-Swift era were due to instrumental biases (e.g. detector 
thresholds effects) rather than to real physical properties of 
the sources. This is the result of a study they performed on a 
sample of 77 Swift GRBs with measured redshifts, where they 
found that those GRBs are on average harder than the pre-Swift 
measurements, i.e. they have higher Ep^ak and lower Ejso- As a 
consequence, they conclude that these correlations cannot be 
used to estimate the GRB redshifts. 

Since no clear conclusion can be drawn from this debate, we 
have also considered this possibility and we calculated all the 
distributions twice: one case where we attribute a pseudo- 
redshift to the 62 bursts which have no redshift, and the other 
case where we randomly attribute one of the 24 secure redshifts 
measured for the HETE-2 long GRBs. Therefore, for all the 
distributions we obtain we will comment on how they change 
between the two cases. 

Throughout the paper we use the following notations: 

- the spectroscopic or photometric redshifts are denoted by z- 

- the pseudo-redshifts are denoted by t 

- Zsampie refers to the value of the redshift used in our study: 
this is the spectroscopic/photometric redshift (z) when it is 
available and the pseudo-redshift (z) otherwise (see Table|2]). 

4. Parameter 1 : Epeak 

4.1. Measuring the Epeak of X-ray flashes 

The peak energy of a burst is generally measured using the 
XSPEC software (Arnaud & Dorman 120031 1 and fitting a phe- 
nomenological model - GRB (Band) Model (Band et al. II 9931 1 
or Cutoff Power Law Model - on the spectral data. However, for 
some X-ray flashes, the peak energy is lower than 10 keV, so 
the data do not constrain the four parameters of the Band func- 
tion. For those bursts, the value obtained is too close to the lower 
boundary of FREGATE (6 keV) and of the WXM (2 keV). In 
this case, we observe a power-law spectrum instead of a Band 
function. To solve this problem, Sakamoto et al. (2004) showed 
that a modified function, called the Constrained Band Model 
( CBM), based on a three-dimensional subspace of the full four- 
dimensional Band function parameter space, could be used. We 
adopted this model to fit the data when it was necessary (see 
Tabled. 

However, in order to check whether such a model could add a 
bias to the estimate of Epeak, we used XSPEC to simulate an 
XRF spectrum as follows: 

- we use the data of a real XRF, in order to have real back- 
ground statistics and instrumental response matrices. 

- we construct a fake burst that follows a Band model with typ- 
ical values of -1.0 and -2.3 for the power-law indices at low 
and high energy respectively. These parameters correspond 
to the mean values for GRB spectra according to BATSE re- 
sults (see e.g. Preece et al. [2000; Kaneko et al. I2006I I. The 
two other parameters, Ep^ak and the normalization, are the 
variable parameters that let us produce fake XRFs of greater 
or lesser intensity. 

- we generate fake data files for XRFs of various brightnesses 
and various Epeak- 

- we then perform a spectral analysis using a Band model and 
a Constrained Band model, and compare the parameters ob- 
tained. 

We find that the Band model is valid for Epeak down to 8 keV. 
For this limiting value, the errors on the spectral parameters are 



slightly larger than those obtained using the Constrained Band 
function, but they are still valid. For Epeak below 8 keV, the 
Constrained Band model is preferable: it gives correct Epeak val- 
ues or provides at least reliable upper limits. 
In all cases we obtain consistent values between the CBM and 
the Band function when both functions are applicable, and these 
values are consistent with the ones we introduced to produce the 
fake XRFs. We thus conclude that the values determined with 
the Constrained Band model are not biased. 

4.2. The observed Epeak distribution 

Several points can be noticed in FigureHJ first, the distribu- 
tion peaks at about 100 keV and extends over the entire energy 
range of HETE-2 (2-400 keV). Second, the two classes of bursts 
(classical GRBs and X-ray flashes) are clearly distinguished. As 
a comparison, the study of Kaneko et al. (2006) on the complete 
spectral catalog of bright BATSE Gamma-Ray Bursts, leads to 
an distribution peaking at about 300 keV, extended at high 
energies up to several MeV and without significant events at very 
low energies. The paucity of soft (~a few keV) and hard (~MeV) 
bursts in our sample is probably not real. For the soft bursts, 
this is clearly visible by a simple comparison of the BATSE 
and HETE-2 Epeak distributions. As BATSE triggers were gen- 
erally done on the count rate between 50 and 300 keV, the ratio 
between the XRFs and the GRBs is low. For HETE-2, the en- 
ergy range extends down to 2 keV and enables the detection of 
more XRFs but still prevents the detection of XRFs with Epeak 
below a few keV. As far as the hard tail of the Epeak distribu- 
tion is concerned, the low sensitivity of the detectors at high 
energy is also one of the causes of the inefficient detection of 
hard bursts. Moreover, two bursts emitting the same energy with 
different hardness do not have the same number of photons: the 
softer GRB has more photons than the harder one. Consequently, 
BATSE-like or HETE2-like missions are probably not able to 
derive a true observed Epeak distribution. The distributions ob- 
tained clearly show that the distribution depends strongly 
on the energy range of the instruments involved. 

4.3. The intrinsic Epeak distribution 

So far, the intrinsic Epeak distribution has been very little dis- 
cussed, mainly due to the difficulty of measuring both the Epeak 
and the redshift of the bursts. On the one hand, BATSE, Beppo- 
SAX and HETE-2 did constrain the £■"'", thanks to their wide 

peak 

energy ranges, but the bursts they detected suffered from a lack 
of redshift determination. On the other hand, thanks to its perfor- 
mance and fast localization capability. Swift currently enables 
more redshift measurements, but the narrow energy range of the 
BAT instrument (15-150 keV) does not allow a good constraint 
of the 

peak 

As mentioned by Amati (120061 1. about 70 intrinsic Epeak are 
available, but this sample contains bursts detected by different 
satellites. Our work is based on a sample containing all the bursts 
localized with a single satellite, hence avoiding selection effects 
inherent to the use of data from different instruments. 
We see in Figure |5]that the fi™^'',;; distribution is broader than the 
El^l^i^ distribution. At high energies, some bursts have an i^™^^ 
reaching a few MeV. Nevertheless, as we previously mentioned, 
the energy range of HETE-2 (2-400 keV) prevents us from draw- 
ing any conclusion about the high energy part of the Epeak distri- 
bution. 
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At low energies, bursts having an ^JJ^^j. lower than 20 keV, and 
defined as X-ray flashes in the observer frame, still have a low 
Epeak in their source frame. For most of those GRBs, the low in- 
trinsic peak energy is due to their intrinsic faintness as they are 
at low or intermediate redshift. (This point is discussed in more 
detail in Section[8]l. 

Joining these two extreme energy ranges, the intrinsic Ep^ak dis- 
tribution of HETE-2 GRBs extends over 3 decades in energy. 
The true width of the distribution is probably even larger be- 
cause bursts with Ep^ak lower than 1 keV are not detected, and 
bursts with Epgak greater than 1 MeV cannot be constrained by 
an HETE2-like mission. 
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Fig. 4. Distribution of the E"**, measured for the complete sam- 

» peak ^ 

pie of 82 GRBs detected and localized by HETE-2. The sub- 
sample's histogram of bursts with secure redshifts (z) is shown 
in red, and the total sample's histogram is shown in grey. The 
bursts on the left part of the histogram correspond to the events 
that have a measured higher limit in of about 4 keV (see 
Table |2]). The energy range of HETE-2 (2-400 keV) is symbol- 
ized by the magenta horizontal solid line. 



4.4. The unbiased Epeak distribution 
AAA . Construction 



Fig. 5. Distribution of for the complete sample of 82 

GRBs detected and localized by HETE-2. The histogram related 
to the sub-sample of bursts with secure redshifts (z) is shown 
in red, and the total sample's histogram is shown in grey. The 
blue dashed histogram corresponds to the studied case 2 (see 
Part |3.3t , i.e. the sample of 62 GRBs without secure redshift that 
are randomly given a redshift from the HETE-2 redshift distri- 
bution. 



the measured E"''^, 

peak 



given m Tableland 100 redshifts ran- 
domly selected within the 90% c.l. error range of their cor- 
responding spectroscopic/photometric redshift or pseudo- 
redshift z,,™. Combined with the 100 simulated . , we 

peak, Sim 

thus obtained a sample of 100 simulated £■'""■, . . 

^ peak, sun 

representing each of the 8200 simulated bursts (b) with a 
normal distribution of the logarithmic values {LE'^^'p: 



Mispj) 



cr 



X exp 



1 (lepj-^LE;- 



(1) 



where le,, ,■ is the variable of the functions fh, Ujm is the 
base 10 logarithm of E'""' 
is set to 0.05. 



, . related to the burst ^ and o", p™. 

peak, Sim p.h 



(2) then, we computed for the 8200 simulated GRBs a weight 
related to their detection significance. For that, we: 



In order to derive the true Epeak distribution, the distance-scale 
is not the only correction that has to be applied to the observed 
Epeak distribution. All the bursts do not have the same brightness 
and are not detected at the same significance in the observer 
frame. We took this into account to correct the intrinsic Epeak 
distribution obtained in the previous section. 
To do this we used the following method: 

(1) in order to take into account the errors on all the pa- 
rameters used (spectral properties and redshift) we have 
computed 'smoothed distributions' by: 

- producing for each burst a set of 100 fi"*^^ randomly se- 
lected within the 90% confidence level (c.l.) error range of 



- determine the maximum redshift Zmax at which the source 
could have been detected by the instruments, by first com- 
paring the SNR computed in Section |2] (see Table |2]i to the 
SNR threshold of the HETE-2 instruments, and then by red- 
shifting the sources until their peak photon flux reaches the 
SNR threshold of the trigger instrument, giving us both the 
visibility distance between the source and the satellite and 
the maximum redshift Zmax- 

- assume that the GRB rate follows the star formation rate. For 
this we have adopted the model SFR2 of Porciani & Madau 



' For some XRFs, the lower error limit could not be constrained. In 
this case the value was set to £°*^j.„„„ = 1 keV. For some other GRBs, it 
was the upper error limit which was not constrained. In this case we set 
the value to = 1 MeV. 

peak^max 
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( 1200 11 1 that reproduces a fast evolution between z = and 
z = 1 and remains constant beyond z>2 



RsfrAz) ^ 0.15/165 



exp{3Az) 
exp(3Az) + 22 



(2) 



derive for each burst the number of GRB per year within its 
visibihty volume 



f 

Jo 



dz 



dV{z) RsfrJz) 



dz 



1 +z 



(3) 



In this equation dV(z)/dz is the comoving element volume, 
described by 



4ndf(z) 



dV(z) _ _c 

dz ~ Tfo (1+ z)2 [Qm(1 + z)^ + Q.k(1 + z)2 + Qa]1/2 



(4) 



where Hq is the Hubble constant, D.k is the curvature contri- 
bution to the present density parameter (Q/f = 1 -Qm - Qa), 
Q.M is the matter density and Q.^ is the vacuum density. 
Throughout this paper we have assumed a flat ACDM uni- 
verse where (Hq, Q.m, Qa)=(65/i65 km s^'Mpc"', 0.3, 0.7). 
This procedure allows us to give each burst a weight (Wb) 
inversely proportional to Nvma.x- The rationale of weighting 
each burst by 1 /Nvmax is the following: the visibility vol- 
ume is different for each GRB of our sample. Moreover, 
each burst observed is randomly taken from all the bursters 
present in its visibility volume. In this way, rare bright bursts, 
having a large visibility volume, will have low weights, 
while faint local GRBs will have heigher weights. This pro- 
cedure also takes into account the fact that the GRB rate 
evolves with redshift, leading to the fact that GRBs are about 
ten times more frequent at z ~ 1 than at present. Note that 
an evolution of the GRB distributions with the redshift is not 
included in this procedure, although such evolution probably 
exists (e.g. Daigne, Rossi & Mochkovitch l2006] l. This case 
is not addressed here because we do not have a sufficient 
number of GRBs. 

(3) the total smoothed distribution (SD) is the sum of the individ- 
ual functions fb(lep,i), normahzed by the corresponding burst's 
weight Wb 



SD 



1 



(LE'°", , ) 



2 Wb X fb(lep,i) 



(5) 



b=l 



with risimu the number of simulations (100) for each burst and 
n bur St the number of bursts contained in our sample (82). 

4.4.2. Results 

The resulting unbiased Epeak distribution for our sample of bursts 
with this method is shown in Figure |6] In this Figure we see 
that the Ep^ait distribution has dramatically changed, with a clear 
domination of the bursts having low intrinsic Ep^ak, i-e. the X- 
ray flashes. 

As explained in Section [3] we tested whether a high accuracy 
for the pseudo-redshifts was necessary, i.e. whether the results 
obtained could depend on these estimated distance scales. To do 
this, we performed the same analysis as described throughout 
this Section, but this time for simulated redshifts randomly se- 
lected within the error ranges based on the initial ones and artifi- 
cially increased by factors of 1.4, 2.0 and 3.0. We find that both 
the distribution of intrinsic Epeak and the unbiased distribution 
of intrinsic Ep^ak do not significantly depend on the accuracy of 





log E;;;;^(keV) 



Fig. 6. Unbiased distributions of the intrinsic Ept,ak for the com- 
plete sample of 82 GRBs detected and localized by HETE-2. 
The smoothed distribution (SD) based on the 20 bursts having 
a spectroscopic/photometric redshift (z) is represented by the 
lower solid red curve and the one corresponding to the com- 
plete sample is shown with the upper solid black curve. The three 
black dotted smoothed distributions show the test performed to 
quantify the impact of the pseudo-redshift accuracy in this study 
(all the pseudo-redshifts are given an error range artificially in- 
creased by factors of 1 .4, 2.0 and 3.0). The blue dashed SD corre- 
sponds to case 2 (see Part |3.3l ). i.e. the sample of 6200 simulated 
GRBs (corresponding to the 62 GRBs of our sample without 
secure redshift) which are randomly given a redshift from the 
HETE-2 redshift distribution. 

the pseudo-z within a factor 2 or 3 (Fig. |6]l. This result can be 
explained by the intrinsic dispersion of the parameters studied 
which is so large that the uncertainties in the distance scales do 
not strongly affect the results. 

Moreover, we also performed this analysis by assuming for the 
distance of each of the 62 GRBs without secure redshift a set 
of 100 redshifts randomly taken from the total redshift distribu- 
tion of HETE-2 long bursts (24 values). The results obtained in 
this case differ significantly because the predominence of low 
E'""'. bursts has disappeared. With this method, the distribution 

peak ' 

is flat and broad. Here again, the distribution's width may not 
reflect the true one, but it contributes interesting information to 
the unbiased Epe^k distribution, which is probably situated - 'in 
reality' - between the two curves presented in Fig. |6] (black solid 
line and blue dashed line). 

5. Parameter 2: Tgo 

In this section we study the duration distribution of the complete 
sample of HETE-2 GRBs. We adopted Tgo as the common defi- 
nition for all the bursts. We derive in the following sub-sections 
the observed, intrinsic and unbiased log Tgo distributions us- 
ing the same method as described for the study of the intrinsic 
log Epeak distribution (Sect.|4]l. 

5. 1 . The observed Tgo distribution 

Tgo is defined as the duration for which 90% of the counts in a 
given energy range are detected. Recall that in this study we only 
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focus on long-duration GRBs. Over the 4 FREGATE bands (A: 
6-40 keV; B: 6-80 keV; C: 30-400 keV; D: >400 keV) we chose 
the one that contained both the classical GRBs and the XRFs, 
i.e. the band B (6-80 keV). The log Tgo distribution in this energy 
range (see Fig.|2ll extends from 0.33 to 2.73. The mean value of 
this distribution is 1.40 and the median value 1.33. 



5.2. The intrinsic T^q distribution 

The histograms are shown in Figure [8] We caution that the ob- 
served Tgo are measured in the same energy range in the observer 
frame (see Table|2|i. Since this energy range depends on the red- 
shift of the source, it is different for each burst and this has an 
impact on the intrinsic duration measurement. 
The histogram of the intrinsic log Tgo in the (6-80 keV) energy 
range extends from -0.26 to 2.40. The mean value of this distri- 
bution is 1.01 and the median value is 1.00. 
Contrary to the log fipeak distribution, the intrinsic log Tgo distri- 
bution is not broader than the observed log Tyo distribution and 
is extended over about the same range (-2.5 decades). 
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Fig. 7. Distributions of the log Tgo measured in the observer 
frame and in band B of FREGATE (6-80 keV) for the complete 
sample of 82 GRBs localized by HETE-2. The histogram related 
to the sub-sample of bursts with secure redshifts (z) is shown in 
red, and that of the total sample is shown in grey. 



5.3. Tlie unbiased Tgo distribution 

Correcting the intrinsic Tgo for the Nvmax (see Part 14.4b we ob- 
tained the unbiased Tgo distribution (Fig.|9]l. 
We note that the GRB population is dominated by events with an 
intrinsic duration of about 10 s (for an energy range of 6-80 keV 
in the observer frame). Nevertheless, a non-negligable number 
of GRBs have a long intrinsic duration lasting about 300 s. 
If we consider the case 2 study (Fig.|9]l, we obtain a significantly 
different unbiased Tgo distribution, more flat than the one ob- 
tained in case 1, i.e. no typical intrinsic duration for the GRB 
appears when considering distances randomly taken within the 
total redshift distiibution of HETE-2 GRBs. 



20 




Fig. 8. Intrinsic log Tgo distributions for the complete sample 
of 82 GRBs detected and localized by HETE-2. The histogram 
of the sub-sample of bursts with secure redshifts (z) is shown in 
red, and the total sample histogram is shown in grey. The blue 
dashed histogram corresponds to the case 2 (see Part l3.3b . i.e. to 
the 62 GRBs without secure redshift that were randomly given a 
redshift from the HETE-2 redshift distribution. 
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Fig. 9. Unbiased distributions of the intrinsic log Tgo for the 
complete sample of 82 GRBs detected and localized by HETE-2. 
The smoothed distribution based on the 20 bursts with spectro- 
scopic/photometric redshifts (z) is represented by the lower solid 
red curve; the one for the complete sample is given by the up- 
per solid black line. The blue dashed SD corresponds to case 2 
(see Part 13.3b , i.e. the sample of 6200 simulated GRBs (corre- 
sponding to the 62 GRBs of our sample without secure redshift) 
which were randomly given a redshift from the HETE-2 redshift 
distribution. 



6. Parameter 3: isotropic energy 

6. 1 . The observed fiuence 

Taking the fluences (i.e. the fluxes integrated over the duration 
of the bursts' spectra) in the HETE-2 Catalog (Vanderspek et 
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al. I2OO8I 1 and applying the same method described in Section |4] 
we show in Figure [10] the histograms of the logarithmic fluence 
in the energy range 2-30 keV. The distribution of the log Sx 
extends between -7.12 and -4.24. The mean and median values 
of this distribution are -5.90 and -6.00, respectively. 

6.2. The isotropic energy 

The equivalent of the fluence measured in the observer frame 
is the energy radiated in the source frame. One easy way to es- 
timate this energy is to assume an isotropic emission (Eiso) in- 
tegrated over a fixed energy-range (1-10"* keV in the observer 
frame, see e.g. Amati et al. 120021 Amati 120061) . We computed 
this quantity for the 82 bursts and obtained the distribution 
shown in Fig.[TT] This distribution of log Eiso extends from 49.50 
to 54.24, and has a mean value of 52.44 and a median value of 
52.62. Two classes of bursts are evident, with a separation at 
about 51.5. Most of the bursts have a log Eiso lying between 52 
and 54 (the classical GRBs) and the second sample has a log Eiso 
lower than 5 1 . 
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Fig. 10. Distributions of the logarithmic fluence in the X energy 
range (2-30 keV) for the complete sample of 82 GRBs localized 
by HETE-2. The histogram for the sub-sample of bursts with se- 
cure redshifts (z) is shown in red, and the total sample histogram 
is shown in grey. 



6.3. The unbiased isotropic energy 

The correction of these Eiso for the Nvmax of each burst gives 
us the unbiased Eiso distribution (see Fig. [T2T l. This distribu- 
tion decreases with Eiso, and the best-fit power-law for the com- 
plete sample distribution, estimated between £,^0=10^" erg and 
£iso=10^"* erg is lO"^''^"* x £',;/^'52^ where £,^0.52 is in units of 
10^2 erg. As for the two previous cases studied, this distribution 
is dramatically different from the observed and intrinsic ones. It 
shows a clear predominence of bursts with low Eiso. This distri- 
bution hence strengthens the result we obtained in Section|4] i.e. 
the predominence of X-ray flashes in the overall GRB popula- 
tion. Here again we can discuss the impact of our second case 
calculation. In this case, the shape of the distribution is similar 
to the one obtained in case 1. 
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Fig. 11. Eiso distributions for the complete sample of 82 GRBs 
localized by HETE-2. The histogram of the sub-sample of bursts 
with secure redshifts (z) is shown in red, and the histogram of 
the total sample is shown in grey. The blue dashed histogram 
corresponds to case 2 (see Part 13.3b . i.e. to the sample of 62 
GRBs without secure redshifts that were randomly given a red- 
shift from the HETE-2 redshift distribution. 



7. Correlations 

Having determined the intrinsic properties of the HETE-2 
GRBs, we test whether any correlations are found between these 
values. We performed a pairwise comparison of the four intrinsic 
parameters {Epgak, Tyo, Eiso, Zsampie) to search for possible corre- 
lations. Six plots are obtained and shown in Figure[T3] The corre- 
sponding correlation coefficients and associated p-values testing 
the hypothesis of no correlation are presented in Table[T] We note 



Table 1. Correlation coefficients (R) obtained for the pairwise 
comparison of the four intrinsic parameters discussed in this pa- 
per {Epeak, T90, Eiso, Zsampie)- The lowcr and upper bounds are at 
a 95% confidence interval. The corresponding p-values test the 
non-correlation by computing the probability of having a corre- 
lation as large as the observed value by chance. 



Tested correlation R 



p-values 



^peak 


0.91« 


1.215x10- 


-32 


<-sample ^ peak 




4.297x10" 


-25 


17 'rinlr 


-0.102+"^ 


3.596x10- 


-1 


'rinir 


-0 354+^-^'^^ 

"-'-'^-0.177 


1.087x10- 


-3 


'T'intr T^intr 
90 ''^peak 
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that the E'^^'J^j-Eiso (Amati) correlation is found (P' panel) and 
that the sub-sample of XRFs fills the gap between XRF 020903 
and the cluster of intrinsic classical GRBs (see next Section). 
This is not really surprising since our pseudo-redshift estimates 
assume that the /s^^^-is/jo relation is valid. It is nevertheless in- 
teresting to note that nearly all of the 82 GRBs in our sample are 
consistent with this correlation with no strong outlier. While it 
is not the purpose of this paper to discuss the reality of E''^'^^^^- 
Zsampie and Eiso-Zsampie Correlations, we note that these correla- 




Fig. 13. Correlations between the intrinsic parameters studied: fi"/^^, Tj^^ ', Eiso and the distance-scale Zsampie- In each panel the 
GRBs that have a measured redshift (z) are represented with the green filled circles, and the ones having a pseudo-redshift (I) 
are symbolized with the blue open circles. In the top left panel, by way of comparison, we show (solid line) the best-fit power-law 
^pmk~'^'^^ Ijo 52 (with in keV and Eiso,52 in units of 10^^ erg) found by Amati (2006 ) and also the vertical logarithmic deviation 
of 0.4 (dotted lines) displayed in Amati (120061 1. In the top right panel the dashed curve marks the redshift evolution of the Ept,ak, 
assuming a burst with e'-^^^i^ - 40 keV (see Stratta et al. |2007l l. 



tions could be partially or totally explained by selection effects 
which practically prevent the detection of low-luminosity GRBs 
at high redshift. In the three last plots, no other 'simple' tight 
correlation is found, but these plots interestingly show the in- 
strumental limits aff'ecting the detection of the bursts (duration, 
energetics and redshift). For instance, we can clearly see that 
HETE-2 was lucky to detect the bright ORB 030329 (see pan- 
els 2 and 6). Hence Figure [13] can be used to show the charac- 
teristics of the bursts missed by HETE-2, i.e. mainly the GRBs 
at low-redshift with high or high Eiso', and conversely the 
GRBs with low Epe^k occuring at high redshift. 



8. The nature of X-ray flashes 

The demonstration of the existence of X-ray flashes (Heise et 
al. 120011 Kippen et al. 1200 It was immediately followed by sev- 
eral theories attempting to explain their nature and to determine 
whether they belong to the same class as the classical GRBs, or 
whether they are different cosmological events. Thanks to vari- 
ous studies of samples of X-ray flashes, several clues now allow 
us to unambiguously associate XRFs with the classical GRBs. 
XRFs are GRBs with lower Epeak and higher fluxes in X-rays 
than in y-rays (e.g. Kippen et al. 120031 12004 . Barraud et al. 
120031 Sakamoto et al. 2005|. The other spectral parameters of 
the prompt emission and the temporal behaviour of XRF after- 
glows at various wavelengths are identical to those of classical 
GRBs (e.g. D'Alessio, Piro & Rossi l2006l l. 
In our complete sample, we have 22 bursts classified as XRFs. 
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Fig. 12. Unbiased Ejso distributions for the complete sample of 
82 GRBs localized by HETE-2. The smoothed distribution based 
on the 20 bursts with a spectroscopic/photometric redshift (z) is 
represented by the lower solid red curve; the one for the com- 
plete sample is shown by the upper solid black line. The blue 
dashed SD corresponds to case 2 (see Part 13. 3t . i.e. the sam- 
ple of 6200 simulated GRBs (corresponding to the 62 GRBs 
of our sample without secure redshifts) which were randomly 
given a redshift from the HETE-2 redshift distribution. The best- 
fit power-law (green dashed line) for the complete sample, esti- 
mated between Eis,, =10^" and 10^"* erg is 10 



-7.94 



(l()"0 



-0.84 



In the following section we discuss some results obtained with 
this sample. 

8.1. The redshift distribution of HETE-2-XRFs 



To explain the origin of the XRFs, Heise et al. (1200 11 1 proposed 
that they are classical GRBs occuring at high redshift. If we con- 
sider the redshift distribution of the HETE-2 complete sample, 
we note that the sub-sample of XRFs has a mean distance scale 
< Zsampie >xRFf^~ 0.82, lower than the mean redshift of the 
'classical' GRBs: < Zsampk >grbs"^= 2.09 (see Fig. [H). The 
mean redshift of the XRFs contained in our sample is also lower 
than the mean redshift of the Swift XRFs derived by Gendre, 
Gain & Piro ( [2007b : < z >xrf'= 1-40 for a sample of 9 XRFs 
with known redshifts. Consequently, we reject this hypothesis, 
as was previously done in other studies, e.g. by Barraud et al. 
( I2003I I who highlighted this fact with the similar duration distri- 
butions they obtained for the XRFs and the long classical GRBs. 



8.2. Do X-ray flashes form a distinct population ? 

We consider in this paragraph the question posed, e.g., in Stratta 
et al. (I2007I I: do the XRFs form a continuum extending the 'clas- 
sical GRBs' to soft energies, or do they form a distinct popula- 
tion, the intrinsic XRFs (i-XRFs)? 

Fig.[T3](panel 1) shows the Amati relation for the complete sam- 
ple of 82 HETE-2 bursts and Fig. [13] (panel 2) plots the 
versus the distance-scales of the sources. We do not find any 
clear evidence for a distinct population and the results indicate 
a continuum between the classical GRBs and the XRFs. Indeed, 



Fig. 14. Four distance-scale cumulative distribution functions 
(cdf). The cdf of the HETE-2 complete sample is shown with the 
green solid line. The two cumulative distributions for the 'clas- 
sical' GRBs and the XRFs are shown with the blue dashed line 
and the red dotted line, respectively. We note that most of the 22 
XRFs are definitely situated at intermediate or low redshift. For 
comparison, we also considered all the long-duration bursts with 
redshifts detected by Swift-BAT from the beginning of the mis- 
sion to GRB 080210 (the last Swift burst with a spectroscopic 
redshift measured to date). The corresponding cdf is the black 
solid line. 



according to the HETE-2 sample, several XRFs are situated in 
the gap between the group of intrinsic classical GRBs and the 
two intrinsic XRFs 020903 (Sakamoto et al. ,2004) and 060218 
(Amati et al. |2007l l, a s is XRF 050416, a soft burst detected by 
Swift (Sakamoto et al. '2006a', '2006b). However, half of the sub- 
sample of the observed XRFs are part of the cluster of intrinsic 
classical GRBs when considered in their source frame. This fact 
was also noted by Stratta et al. (I2007I I. 

9. The rate of GRBs in the Local Universe 

9.1. The rate of GRBs in the Local Universe, as measured by 
HETE-2 

The various corrections (distance scale, number of GRB within 
the visibility volume) applied to the distributions obtained in the 
observer frame allowed us to obtain the unbiased distributions 
of the intrinsic parameters, i.e. distributions that reflect more 
closely the true distributions of these parameters and the dom- 
inant characteristics of the gamma-ray burst population. 
This study also allows us to derive the rate of GRBs detected and 
locaUzed by HETE-2 (R^^). For that purpose, we consider that 
each GRB in our sample contributes to the local rate in propor- 
tion to: 



hi, - 



■^Vmaxiz — Zb,max) 



(6) 



with the Nviz) computed according to Eq. [3] and we obtain the 
rate of HETE-2 GRBs during the mission. 



y,oc 



"Inirsi 



(7) 
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which is 3.30 Gpc ^. 

In order to normaUze this rate per year, we took into ac- 
count the effective monitoring time of the WXM, obtained 
from: 

rj. _ 7"mission X 

^ m — ^ '-> cov \o) 

An 

where T,,, is the effective monitoring time of the WXM, 
T mission - 69 months is the duration of the mission, =37% 
is the mean observation efficiency during T mission, and S cm- - 
2n{\ - cos 45°) = 1.84 sr the sky-coverage of the WXIvfl 
The effective monitoring time of the WXM is hence 
Tm — 0.31 yr Using this, the rate of GRBs in the Local Universe 
per Gpc^ and per year can be found to be ~10.6 Gpc^^yr ' . This 
result is a lower limit because as shown in Section H) the bursts 
that dominate the overall population are the ones with low 
and as we showed in Section |7] HETE-2 missed the bursts with 
intrinsic Epeak lower than 1 or 2 keV as well as bursts occuring 
at very high redshifts. 

If we consider the unbiased distributions obtained in case 2 
(shown as dotted blue lines throughout the different Figures), we 
obtain 1.02 Gpc"^, corresponding to a rate of ~3.3 Gpc^^'yr"'. 
As shown in Section H] the population of gamma-ray bursts 
is dominated by the X-ray flashes. This is understandable as 
XRFs are soft but also faint in the observer frame, according to 
the hardness-intensity relation derived by Barraud et al. (2003). 
Therefore, if the rate of detected bursts is in fact higher for clas- 
sical GRBs than for XRFs, we can guess that this is only due to 
instrumental limitations, such as the size and the sensitivity of 
the detectors. 



and 1.1 Gpc yr . This rate is in agreement with the one found 
by Schmidt (20011: 0.5 Gpc^^yr ' with 13 91 GR Bs contained 
in the BATSE DISCLA 2 Cat alog (S chmidt |T999] l and the one 
recendy found by Liang et al. (|2557] | (1.12+°^3 Gpc^^yr ') with 
a sample of 45 HL Swift-GRBs. The rate of HL HETE-2 GRBs 
is also consistent with all these results, since we obtain a lower 
limit Rf- >0.45 Gpc ^r ' (see PartlQ. 

(2) GRB 060218, detected by Swift-BAT (Cusumano et al. 
,2006) , was found to be at very low-redshift z =0.0331 (Mirabal 
et al. I2006al l. Moreover, its association with the supernova 
SN 2006aj (Pian et al. '2006) recalled the characteristics of 
GRB 980425, the closest GRB ever identified, at a redshift 
z =0.0085 (Tinney et al. [T998] l. These GRBs associated with 
SNe are peculiar in the sense that they constitute the two nearest 
long duration bursts detected in eight years. Liang et al. (2007]l 
show that the detection of these two events implies a local rate 
of low-luminosity GRBs (LL-GRBs) of Kll ~ 800 Gpc^^yr 
in agreement with Guetta & Delia Valle ( I2007I I who find a rate 
of GRB 980425-like events of ~ 430 Gpc ^r"' and a rate of 
GRB 060218-like events of ~ 350 Gpc-^yr '. 
The rate of LL-GRBs, which is considerably higher than the 
rate of HL-GRBs, could be explained in two ways: first, the 
log(N) - log(P) disti-ibution of BATSE bursts often used to fit 
the GRB luminosity function cannot constrain it down to such 
low luminosities, because a GRB 060218-like event could not 
trigger BATSE. Second, the use of a phenomenological model 
for the GRB LF requires a choice for the extrema values Lm/„ 
and L„,ax- However, due to instrumental limitations, the true Lm,„ 
of GRBs is unknown. Thus, the number of low-luminosity GRB 
is quite uncertain. 



9.2. Comparison with tine rates of GRBs obtained in previous 
studies 

Several studies have derived local rates of GRBs using various 
methods mainly based on theoretical considerations. 
According to the results obtained so far, the predicted rate of 
GRBs in the Universe seems to be between unity (or less) GRB 
and several hundreds per Gpc^ and per year. These results seem 
to be related to the nature of the bursts considered in the studies, 
i.e (1) 'classical GRBs' also called 'high-luminosity GRBs' 
(HL-GRBs) and/or (2) sub-luminous (low-luminosity) GRBs 
(LL-GRBs). 

(1) Several authors obtained a GRB rate close to 1 Gpc"^ yr 
The method they used is generally based on the determination 
of the luminosity function (LF) by fitting a phenomenological 
model to the log{N) - log(P) relation of well-known GRB 
catalogs. Assuming that GRBs trace the star formation rate 
(SFR), a model of the SFR is adopted for it and the rate of GRBs 
is determined. 

For instance, Guetta et al. (|2004l |2005] l use the 2204 GRBs of 
the GUSBAD Catalog (Schmidt '2004) and a sample consisting 
of 595 long GRBs detected by BATSE. Following Schmidt 
( I2001I I, they use these samples to find the luminosity function 
by fitting a single power-law model to the BATSE observed 
peak flux distribution. In addition, they assumed that the 
GRB rate follows the SFR and hence adopted either the SFR 
model SFR2 of Po rciani & Madau (2001) or the SFR model of 
Rowan-Robinson ( 119991 1. The estimated rate lies between 0.1 

Recall that throughout this study, we only consider the GRBs local- 
ized by the WXM. 



9.3. Discussion 

We showed that the study of the complete sample of GRBs 
localized by the WXM on-board HETE-2 led to a rate 
Rq^ <: 10.6 Gpc^^'yr Our result is independent of any assump- 
tion on the luminosity function, as it is based on measured data. 
We note that this measured rate is at least 10 times higher than 
the HL-GRB rate and 10 times lower than the LL-GRB one. How 
can we explain this 'intermediate' result? 

- (R^^ < Rq^): We first note that we have no event like 
GRB 980425 or GRB 060218 in our sample. Previous 
studies showed that it was impossible for an HETE-like 
mission with a WXM-like camera to detect a burst like 
XRF 060218. This burst occured at a distance d =145 Mpc, 
and Soderberg et al. (2006d) show that the WXM could 
have detected it at a maximum distance of 110 Mpc. For a 
burst similar to GRB 980425 occurring at li = 36.1 Mpc, 
the WXM could have detected it to a distance d - 60 Mpc. 
However, Guetta et al. (2004) show that with their model of 
the LF, the local rate of nearby bright bursts observable by 
HETE-2 within this distance is only 0.057 yr 



- (R^ > Rq ): If we now imagine that HETE-2 was a mission 
that could not detect X-ray flashes, we can try to re-do all the 
steps to derive the rate of GRBs detected and localized by it 
(see Part. 19. II ) in the sample of 60 'classical' GRBs. We ob- 
tain R^- <: 0.45 Gpc^^yr"'. This result is in good agreement 
with the rates usually obtained when only the HL-GRBS are 
considered: ~[0.1-1] Gpc^^'yr"'. 
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Consequently, we infer that the difference between the rate 
of HL-GRBs obtained with BATSE and the rate measured 
with HETE-2 is due to the X-ray flashes. These soft events 
cannot be detected in the BATSE energy range (50-300 keV). 
With our method based on observed properties of GRBs, 
XRFs and classical GRBs can be distinguished. The average 
lower redshift of XRFs leads to higher weights for the XRFs 
than for classical GRB^. In fact, it appears that XRFs con- 
stitute an intermediate class in terms of frequency between 
the HL-GRBs and the LL-GRBs, and the rate of XRFs could 
be higher than the rate of classical GRBs by a factor of ~10-100. 

Despite their higher softness, faintness and frequency, our 
study brings the XRFs closer to the classical GRBs. This is not 
the case for GRB 980425 or XRF 060218-Uke events that are 
undetectable and out of reach of a HETE-2-like mission. 

9.4. Comparison with fhe rate of Type Ib/c SNe 

If we compare the measured rate of GRBs based on the com- 
plete HETE-2 sample and the rate of SNe Ib/c considered in 
Soderberg et al. ( I2006db . 7?^^"*^" = 9+ ^ x 10 ^ Gpc-^yr ' mea- 
sured by Cappellaro, Evans & Turatto ( 1 19991 1 and Dahlen et al. 
( |2004] |. the ratio is >0.1%. 

This low ratio highlights the fact that not all the type Ib/c SNe 
produce GRBs. The main explanation invokes a mildly relativis- 
tic jet, which is the source of the GRB, but that does not appear 
in all SNe Ib/c. Indeed, if the optical luminosities of GRBs and 
SNe can be considered to be similar, optical spectroscopy reveals 
broad absorption lines in SNe spectra associated with GRBs, in- 
dicative of fast ejecta. These peculiar highly energetic SNe, often 
designed Hypernovae (HNe) would not represent more than 3% 
of the SNe Ib/c, as obtained by the studies based on the compar- 
ison of radio luminosities of GRB afterglows and SNe Ib/c (see 
e.g. Berger et al. |2003t Soderberg et al. |2006bt Soderberg [20U7l ). 
From previous studies of the predicted local rate of GRBs, and 
with our measured local rate, we now have a ratio between the 
SNe Ib/c and the GRBs lying in the range ~[0.1%-3%]. This 
small ratio could imply that the escape of the ultra-relativistic 
jet from the progenitor star is not the only ingredient needed to 
produce a GRB. In addition, several studies of stellar rotation, 
binarity, asymmetry and metallicity have shown that the progen- 
itors must have some special characteristics in order to produce 
a gamma-ray burst. 

10. Summary 

In this paper we have presented the study of a sample con- 
taining all the gamma-ray bursts detected by HETE-2 above a 
certain threshold and localized by the WXM. Our main goal 
was (1) to provide intrinsic distributions of the main global 
properties associated with the GRBs, i.e. taking into account 
the cosmological effects. The parameters studied are the peak 
energy {Epeak), the duration (Tgo) and the total energy {Eiso). 
The second purpose of this paper is (2) to derive unbiased 
distributions of these properties, taking into account the spatial 
density of GRBs, and to measure the rate of HETE-2 GRBs in 
the local universe. 

(1) The main conclusions are the following: 

' Recall that the weight taken for each burst is inversely proportional 
to the number of GRBs within its visibility volume and is obviously 
redshift-dependent. 



- the intrinsic Epeak distribution is quite broad, containing at 
the extreme sides of the energy range a large sample of X- 
ray flashes and few 'classical' hard bursts. Nevertheless, we 
caution that the true width of this distribution is probably 
even broader than derived here, since the detector character- 
istics prevented HETE-2 from detecting bursts with an in- 
trinsic peak energy of 1 or 2 keV or, at the opposite end, 
bursts with an intrinsic peak energy of a few MeV. 

- the intrinsic Tgo distribution shows that HETE-2 did not lo- 
calize very long duration bursts (Tgo >1000 s). This may 
be a limitation of the orbit because due to its antisolar posi- 
tion, the detectors are turned off every half-orbit when they 
have the Earth in their field of view, i.e. every ~ 45 min. On 
the other hand, some events have an intrinsic duration below 
2 seconds and might be considered as 'intrinsically short' . 

- the intrinsic Ejso distribution extends over more than four 
decades in energy. Most of the observed bursts have Eiso 
higher than 10^^ erg: the 'classical' GRBs. The distribution 
extends down to lower Eiso (10^^-10^' erg), corresponding 
to the sample of HETE-2 XRFs. 

(2) The main conclusions are the following: 

- the 'unbiased' distributions are dramatically different from 
their uncorrected shapes. Indeed, the unbiased Epeak distri- 
bution shows a predominence of bursts with intrinsic Ep^ak 
lower than 10 keV and the unbiased Eis,, distribution shows 
a predominence of bursts with Ei^,, lower than 10^"- 10^' erg. 
The Eiso distribution can be fitted between 10^" and lO^'* erg 
by a simple power-law with a slope of -0.84. 

However, we caution that the 'simple' distributions of the 
intrinsic parameters probably suffer from biases. The same 
point may be made for the shapes of the 'unbiased' distri- 
butions: the true ones are probably also different, as shown 
by the two studies performed using different distances for the 
bursts without secure redshifts - case 1 : pseudo-redshifts and 
case 2: known spectroscopic redshifts randomly attributed 
- which lead to different shapes for the Epeak and the Tgo 
distributions. This emphasizes the predominence of X-ray 
flashes and the existence of a typical intrinsic duration for 
the GRBs in the first case, whereas a nearly constant number 
of GRBs through both the Epeak energy range and the du- 
rations is highlighted in the second case. We may guess that 
the true distributions are situated between the results of these 
two cases. 

- a possible application of our study is the determination of the 
rate of GRBs measured by HETE-2 in the Local Universe. 
We obtained a lower limit of 10.6 Gpc^^yr \ which is in- 
termediate in magnitude between the rates usually found 
when a population of High-Luminosity GRBs is considered 
(~1 Gpc^^yr ') or when the Low-Luminosity GRBs are also 
taken into account (^100 Gpc^^yr '). We explain this result 
with two arguments. First, HETE-2 did not detect LL-GRBs 
like GRB 980425 or XRF 060218 - mainly due to the small 
size of its detectors - so the rate cannot be as high as the one 
found for LL-GRBs. Second, it is the sample of XRFs de- 
tected by HETE-2 that makes the difference between the rate 
we obtained and the one found for HL-GRBs. As a conse- 
quence, thanks to the construction of unbiased distributions 
of the intrinsic properties of HETE-2 GRBS, we experimen- 
tally showed in this paper the dominance in terms of fre- 
quency of the X-ray flashes over the 'classical GRBs', a fre- 
quency which is nevertheless lower than the sub-luminous 
bursts, and much smaller than the Type Ib/c supernovae. 
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Table 2. Observed properties and distances used in this study. 

The GRB names in the format YYMMDD are given in column 1. Columns 2 to 5 contain the spectral parameters measured for each burst: the 
slope of the power-laws at low-energy (a) and at high-energy when possible (/?); the peak energy (fi"*^^-) in keV and the fluence in the energy 
range 2-30 keV (Sx), in units of 10"' ergs cm"'. The Tgo durations in the FREGATE band B (6-80 keV) are reported in column 6. The signal- 
to-noise ratio (SNR) corresponding to the best combination of energy range and time resolution (see Section|2} are in column 7. The Zsampie (i-e. 
the spectroscopic/photometric redshifts when available or the pseudo-redshifts otherwise) are reported in column 8. The last column contains the 
references for the redshifts. 
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581.25 
18.02 

9.26 
13.25 
193.85 
28.11 

9.74 
88.56 
21.43 
15.49 
12.84 
14.23 

8.46 
78.68 

8.7 
151.81 

5.81 
92.35 



< 



0.17 

034+0.22 

O lQO+(i'.^^4 

-0.004 
1 Q+0.2 

-0^4 

1 98+*^-' 

^•^°-0.7 
2_n+0.5'5 



2.21 



0.7 
73+°-i3 

1.21+1 

3.48+|-^ 
1 9CC+0.005 

i.Z,JJ_„ 

41+°-°' 

-fi-fii 

9 T 1 +0.01 

^-'^-o.pi 

2-8!!| 
1 +0.5 
^-0.5 

1 l+o^ 

4 15+1-** 
1 006+°-"°i 

2 2+ 

1 98+'''-^2 

^■^°-0.02 

> 2.3 

1 5216+"-™'' 
n ifi8s+ooooi 

ll+°-°4 

"■^^-0.04 

3.07+!-' 

2.6584t 
0.86i! 

n 787+0 001 

54+°- 15 

"-'^-0.15 

0.89!°-2 

1 77+0:5 

' -0.5 
0.83+°-3 



0.26+ 



6.04 
6.67 

94+*^:^ 
2 05+°-** 
2 14+°-'* 

ii+°-i' 

2 17+°-' 

^•^ -0.7 



[0] 
[0] 
[0] 
[0] 
[0] 
[0] 
[0] 

[i;2] 

[0] 
[0] 
[0] 
[0] 
[3] 
[4] 
[0] 
[0] 
[0] 
[0] 
[0] 
[0] 
[5] 
[6] 

[7; 8; 9] 

[0] 

[0] 

[0] 

[0] 

[10; 11] 
[12] 

[13; 14; 15; 16; 17] 
[0] 

[18; 19;20;21] 

[22; 23; 24] 

[0] 

[0] 

[25] 

[0] 

[26] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 
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Table 2. (continued) 



GRB 


a 


P 


(kev) 


Sx 


T90.B 


SNR 


Zsample 


■^^f- ^sample 



-6.75 



031220 

040319 

040423 

040425 

040511 

040701 

040709 

040802 

040825A 

040825B 

0409 12A 

0409 12B 

040916 

040924 

041004 

041006 

041016 

041127 

041211 

050123 

050209 

050408 

050509 

050729 

050807 

050922 

051021 

051022 

051028 

051211 

060115 

060121 



-i.ji_02 

-0.89!»:f 
-1.02+!}3i 

-0.86: 
-0.92 



-3.25 



-0.14 



_i 9 c -1^0.08 

-0.85+»f 
-1 51 +"-22 
-1 48+''-2i 

-1 17+'-54 
^ ■ ^ ' -0.62 

_i 9C-1-O.82 
^■^-"-0.58 
-1 

-1 03+""'^ 
^■"-'-'-0.09 

-1 3+11^ 

_i T-fO.05 
^■■-'-0.05 
1 1 T -1^0.36 
^■^■-'-0.21 



-0.21 

-1 01 +"-27 
^■^^-0.16 

-0 69+''-4^ 

u. /i _() 2 

-1.52+"-^ 

-1 76+^-'''^ 
^ ■ ' "-0.O6 



_0 61+0-43 
^■"^-0.32 
1 <:t -1-0.42 



-0.22 

_ 1.04+"- 11 
-1 01 +^■'"'2 

^■^^-0.03 

_0 94+0.18 

^■^^-0.16 

-0 07+"-^ 

^■"-"-0.41 

-0 79+"- '2 
'^•'^-o.ii 



9 T-i-0.25 
"^■■-'-0.31 



_o 49+0.24 
^■^^-0.4 



-2.3 

_9 97-1-O.24 
^■^'-0.47 

-1 48+0-16 

^■^"-0.54 



-2.3 

_9 9-1^0.16 
^■^-0.28 



-2.3 
-2.3 



1 QC-^0.25 
^■^-'-0.14 



1 81+012 



56.6+1^ 
30.7+f 

9QQ Q+72.9 
^^^■^-48.7 
Qo Q+15.7 
^■^■^-11.3 

< 3.44 
70.9+^-3 
92.2+jf 

60!^« 
25.1!^^/ 

^ '-13 

< 3.5 

41 1+2-3 

53-7!"/ 
47.7!i, 
165.3!^3i. 

35+"' 
■^-'-10.2 

1 09 + 1 16 
^-''^-76.3 

40.2+^ 
445+555 

^^■^-350.5 

25.9+i''35 

< 19 

69.47+J82? 
130.5 

213!!^ 

94.Q+167 
121+33 

94.6; 



-31.12 
c-i-50.9 
-"-26.8 
fS+26 . 



-20.3 
4-21.4 
-17.6 
'•^.2 



5 5+"-'' 
•'-'-0.5 

"-0.5 

99 7-F2.6 

-2 2 

17.2+j:^ 
31.7+"-3 

5 44+(f55 
-■■^^-0,55 

80.6+f 

4 7Q-f0.68 
^- ' ^-0.68 

9 6+' -^ 

1 9+0.9 
^^-0.9 

20.36+2"4 

9 5+"-^^ 
^■•'-0.95 
7 74+0.91 
'■'^-0.91 
90 4+0.24 
^-"■^-0.24 

+ 19.7 
-18.4 
+0.6 
-0.9 
H.53 
-1.93 

26.9!?"^ 

10 6+" '* 

"■'-3.7 

-r ^+0.82 

-'■-'-5.4 

32.4+'^-3 

g+0.6- 

-0.6 
T Q7+0.69 

■^■^ -0.65 

14 15 + '*^'* 

c 4+0.54 

-'■^-0.54 

214+2 
7 QQ+0.59 

+0228 
0.765 n228 
99 4+0.4 
^^■^-0.5 

^■98tg 



38 9+"-* 

8.21_,no 



9.69 + 0.93 
6.1+0.85 
45.87 ± 1.13 
138.81 + 1.65 

45.8 ±0.53 

11.67 + 5.72 
82.48 ± 0.24 
3.05+0.18 
39.2 + 2.01 
15.8 + 5.45 
9.21+0.27 

122.41+7.18 
349 + 10.23 
3.37 + 0.08 

50.13 + 2.49 
22.08 + 0.33 
21.96+ 1.7 

49.14 + 2.26 
113.84 + 4.2 
14.72+ 1.17 

18.38 + 0.59 

28.39 + 0.56 
21+5 

5.23+0.51 
10.32+1 
6.13+ 1.4 
36.79+ 1.13 
178 + 8 
14.66+1 
4.9 + 0.61 

20.68 + 0.45 
2.61+0.1 



13.9 
23.12 
17.47 
11.16 
36.61 
12.76 
58.3 
39.39 
10.47 
10.3 
39.44 
8.19 
6.66 
85.97 
35.64 
120.37 
13.59 
25.19 
44.38 
25.06 
15.37 
50.85 
7.06 
24.35 
19.17 
38.25 
35.54 
226.38 
25.35 
18.83 
87.09 
80.37 



i.jj j 2 

1 79+''2 

1.26+«-^ 

2.23^^ 

1.83!«:^ 

0.2146+:- 
1 +0.3 
^-0.3 

1 01 +0.65 
1-81-0.65 

1 06+"-^ 

9 9+O.S 
^■^-0.8 

33+°-3 

1 563+°"°^ 
< 0.7 

859+° °°^ 

'-'■°-'^-0.005 

0-58+JJ-2 

Q7jg+0.005 
^■'^"-0.005 

3.49+j-^ 
0.96+"! 
3 29+0-^ 

1 co+0.85 
^■-^■'-0.85 
9 Q^ + 1.6 

1 2357+"' 
0.68+;}3 

9 c+0.75 

-0J5 

79+'^ " 
^■'^-0.75 

9 1Q0+0.001 
^■^^°-o.ooi 
1.37+°-^ 
n o+o.?if 

'-'■'5-0.01 

3.66+|-^ 

4.83!;| 
^52+0.25 

1 92+^-^^ 



.,.j002 
0.0002 
-0.3 
^-0.3 



[0] 

[0] 

[0] 

[0] 

[0] 

[27] 

[0] 

[0] 

[0] 

[0] 

[0] 

[28] 

[0] 

[29] 

[0] 

[30] 

[0] 

[0] 

[0] 

[0] 

[0] 

[31; 32] 

[0] 

[0] 

[0] 

[33] 

[0] 

[34] 

[0] 

[0] 

[0] 

[0] 



[0]= our work; [1]= Djorgovski ( I2OOII ; [2]= Price et al. ( l2002l l: [3]= Hjorth et al. j2003bb : [4]= Berger et al. ( l2007t : [5]= Barth et al. i2003] |; [6]= 
Jakobsson et al. ( I2005ab : [7]= Savaglio et al. ( [2002t : [8]= Castro-Tirado et al. ((2002); [9]= Lazzati et al. f2006 l): [10]= Vreeswijk et al. ^2M2\ : 
[11]= Fox et al. (|2003); [12]= Levan et al. l [2006t : [13]= Ando et al. (2003ali; [14]= Ando et al. (2003b}; [15]= Greiner et al. (2003a); [16]= Price 
et al. ( |2003t ; [17]= Shin et al. (2006); [18]= Martini et al. ( |2003] ); [19]= Rol et al. ( |2003l ; [20]= Fugazza et al. (2003); [21]= Maiorano et al. 
( I2006t ; [22]= Greiner et al. (|2003bj); [23]= Caldwell et al. ( |2003] |; [24]= Bloom et al. 12003); [25]= Jakobsson et al. ( .2004 J ; [26]= Rau, Salvato 
& Greiner f2005) ; [2 7]= K elson et al. ( l2054t ; [28] = Stratt a et al. ( 12507} ; [29]= Wi ersema e t al. ( |20M] |; [30]= Fug azza e t al. ( 12504} ; [31]= Berger, 
Gladders & Oemler ( 12505} ; [32]= Prochaska et al. f2055} ; [33] = Jakobsson et al. ( |2005b| l; [34]= Gal- Yam et al. ( |2505} . 
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Table 3. Intrinsic properties of our complete sample of HETE-2 localized GRBs. 

For each burst, the values reported in columns 2, 3 and 4 are base 10 logarithms which enable us to obtain the individual lognormal distributions 
of the Epgak< T9o,B and The z^ax used to derive the corresponding V^,^ are shown in column 5. The resulting V/V,^^ are reported in the last 
colurrm. 

~GRB log^pl't iQg^^QM log£»o.52 Zmox V/Vma. 

001226 1.648 + 0.091 1.66 + 0.037 0.516 + 0.112 1.12 + 0.19 0.779 + 0.012 

010213 0.593 + 0.03 1.239 + 0.028 -2.172 + 0.023 0.39 + 0.08 0.091 +0.009 

010225 2.003 ±0.122 0.403 ± 0.059 -0.173 ±0.08 1.5 ±0.3 1 

010326 2.341 ±0.111 -0.264 ±0.106 0.518 ±0.115 4.85 ± 1.47 0.552 ±0.049 

010612 3.306 + 0.168 0.916 + 0.095 1.457 + 0.11 9.22 + 2.2 0.61 +0.052 

010613 1.891 +0.087 1.952 + 0.034 0.711 +0.12 1.42 + 0.28 0.173 + 0.016 

010629 1.938 + 0.06 0.925 + 0.053 0.357 + 0.054 1.86 + 0.49 0.228 + 0.029 

010921 2.128 ±0.053 1.177 ±0.014 -0.076 ±0.028 1.18 ±0.01 0.104 ±0.001 

010928 3.293 ±0.073 0.841 ±0.048 1.814 ±0.12 9.23 ± 1.43 0.406 ±0.027 

011019 1.525 + 0.178 1.168 + 0.089 -1.324 + 0.123 0.65 + 0.34 0.63 + 0.038 

011130 0.43 ±0.162 0.958 ±0.068 -2.26 ±0.06 0.19 ±0.05 0.705 ±0.038 

011212 0.471 ±0.177 1.476 ±0.063 -1.703 ±0.061 0.38 ±0.16 0.679 ±0.028 

020124 2.531 ±0.053 1.086 ±0.008 1.999 ±0.083 5.58 ±0.06 0.532 ±0.005 

020127 2.481 ±0.082 0.399 ±0.032 0.553 ±0.059 3.45 ±0.38 0.376 ±0.016 

020305 2.995 + 0.17 1.118 + 0.058 1.22 + 0.107 4.11 +0.91 0.367 + 0.03 

020317 1.98 ±0.093 0.36 ±0.06 0.23 ±0.185 2.59 ±0.41 0.727 ±0.016 

020331 2.48 ±0.077 1.755 ±0.057 1.189 ±0.073 3.83 ±0.74 0.455 ±0.032 

020625 1.141 +0.135 0.991 +0.096 -1.498 + 0.126 0.53 + 0.29 0.671 +0.035 

020801 2.054 + 0.129 2.339 + 0.117 1.241 +0.222 1.68 + 0.85 0.507 + 0.049 

020812 2.752 + 0.182 0.606 + 0.103 1.063 + 0.121 4.78 + 1.39 0.716 + 0.035 

020813 2.773 ±0.042 1.588 ±0.002 1.982 ±0.016 3.8 ± 0.03 0.169 ±0.002 

020819 1.858 ±0.079 1.304 ±0.067 -0.251 ±0.076 0.96 ±0.02 0.124 ±0.002 

021004 2.492 ±0.104 1.17 + 0.013 0.754 + 0.087 3.1 +0.02 0.655 + 0.006 

021016 3.037 + 0.17 1.358 + 0.115 1.517 + 0.096 5.01 + 1.83 0.451 +0.064 

021021 1.56 + 0.16 0.961 +0.102 -0.751 +0.143 1.08 +0.3 0.927 + 0.009 

021104 1.84 ±0.109 1.002 ±0.06 0.097 ±0.116 1.25 ±0.21 0.77 ±0.012 

021112 2.517 ±0.148 -0.216 ±0.137 0.482 ±0.168 4.74 ± 1.22 0.875 ±0.015 

021211 1.977 + 0.03 0.321 +0.015 0.348 + 0.043 2.93 + 0.02 0.149 + 0.002 

030115 2.484 + 0.094 0.802 + 0.043 0.592 + 0.03 4.04 + 0.11 0.436 + 0.006 

030226 2.481 +0.053 1.409 + 0.014 0.923 + 0.065 2.9 + 0.03 0.556 + 0.004 

030324 3.237 ±0.255 0.374 ±0.097 0.847 ±0.035 6.05 ± 1.26 0.586 ±0.032 

030328 2.515 ±0.029 1.739 ±0.005 1.381 ±0.013 4.08 ±0.02 0.226 ±0.001 

030329 1.915 + 0.008 1.345 + 0.004 0.03 +0.008 1.16 + 0.00 0.007 + 0.000 

030416 0.295 ±0.173 1.1 ±0.042 -2.501 ±0.139 0.21 ±0.03 0.228 ±0.01 

030418 2.309 ±0.147 1.558 ±0.113 1.175 ±0.055 3.19 ±1.07 0.906 ±0.015 

030429 2.122 ±0.069 0.554 ±0.054 0.54 ±0.049 3.39 ±0.02 0.718 ±0.004 

030519 2.409 ±0.036 0.841 ±0.018 1.21 ±0.067 3.74 ±0.29 0.075 ±0.005 

030528 1.753 + 0.043 1.549 + 0.02 0.483 + 0.035 1.34 + 0.01 0.307 + 0.003 

030723 0.793 ±0.279 0.795 ±0.045 -1.258 ±0.105 0.6 ±0.09 0.802 ±0.013 

030725 2.282 ±0.045 1.965 ±0.036 0.752 ±0.016 2.6 ±0.33 0.132 ±0.012 

030821 2.372 + 0.063 0.845 + 0.048 0.719 + 0.058 3 +0.53 0.447 ±0.024 

030823 1.717 + 0.071 1.439 + 0.043 0.101 +0.054 1.14 + 0.23 0.53 + 0.022 

030824 0.74 + 0.17 0.915 + 0.041 -1.659 + 0.089 0.29±0.1 0.551 ±0.03 

030913 2.996 ±0.147 -0.012 ±0.137 1.023 ±0.189 8.19 ±2.21 0.772 ±0.032 

031026 3.595 ±0.234 0.947 ±0.093 2.178 ±0.136 6.89 ±1.69 0.969 ±0.004 

031109A 2.564 + 0.044 1.47 + 0.024 1.02 ±0.025 2.93 + 0.38 0.129 + 0.011 

031109B 2.231 +0.174 2.145 + 0.129 0.664 + 0.099 3.09+ 1.18 0.94 + 0.01 

031111A 3.116 + 0.045 0.402 + 0.037 1.288 + 0.053 10.51 + 1.17 0.182 + 0.018 

031111B 0.752 ±0.272 1.359 ±0.027 -1.808 ±0.14 0.19 ±0.05 1 

031203 2.683 ±0.074 0.517 ±0.057 1.661 ±0.057 6.63 ± 1.31 0.26 ±0.034 

031220 2.204 + 0.19 0.585 + 0.117 0.013 + 0.062 2.08 + 0.82 0.64 + 0.048 

040319 2.183 + 0.129 0.37 + 0.124 0.5 +0.13 2.9+ 1.31 0.421 +0.062 

040423 1.845 + 0.085 1.314 + 0.078 0.607 + 0.102 1.8 + 0.58 0.507 + 0.041 

040425 3.001 ±0.064 1.634 ±0.038 1.307 ±0.07 2.73 ±0.35 0.739 ±0.009 

040511 2.421 ±0.062 1.22 ±0.05 1.255 ±0.076 3.46 ±0.63 0.366 ±0.03 

040701 0.405 + 0.146 0.955 + 0.122 -2.096 + 0.074 0.3 + 0.01 0.384 + 0.013 

040709 2.15 + 0.051 1.618 + 0.041 0.898 + 0.035 2.49 + 0.5 0.182 + 0.021 

040802 2.43 +0.077 0.033 +0.062 0.955 +0.115 3.74 + 0.79 0.358 + 0.036 
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GRB loggg^-, log 7-^03 logg/.».52 z„a. VrV, 



040825A 


2 056 


+ 


0.133 


1.287 + 0.068 


-0.088 + 0.069 


1.17 + 0.35 


0.772 + 0.02 


040825B 


1 942 


+ 


0.123 


0.701 + 0.098 


0.698 + 0.039 


2.49 + 0.52 


0.838 


+ 0.014 


0409 12A 


1 247 


+ 


0.095 


0.851 +0.064 


-0.782 + 0.24 


0.61 +0.36 


0.162 


+ 0.043 


0409 12B 


1 62 




0.194 


1.68 ±0.015 


0.327 ±0.131 


1.82 ±0.01 


0.754 


± 0.007 


040916 


0.465 


+ 


0.163 


2.401 ± 0.052 


-1.329 ± 0.052 


0.42 ± 0.18 


0.856 


± 0.028 


040924 


1.88 : 




0.014 


0.259 + 0.006 


0.292 + 0.029 


2.5 + 0.02 


0.129 


+ 0.001 


041004 


1.938 


+ 


0.059 


1.504 + 0.037 


0.632 ± 0.235 


1.13 + 0.25 


0.213 


+ 0.019 


041006 


1.915 


± 


0.016 


1.11 ±0.004 


0.219 ± 0.013 


2.37 ± 0.02 


0.092 


± 0.001 


041016 


3.263 


± 


0.288 


0.689 ± 0.096 


1.64 ±0.243 


4.95 ± 1.35 


0.7 ±0.031 


041127 


1.809 


+ 


0.091 


1.417 ±0.067 


0.458 ± 0.08 


1.53 ±0.49 


0.341 


± 0.038 


041211 


2.804 


+ 


0.188 


1.42 + 0.054 


1 .496 + 0.263 


8.66 + 1.5 


0.389 


+ 0.029 


050123 


2.013 


± 


0.092 


0.765 ± 0.085 


0.434 ± 0.142 


2.73 ± 0.86 


0.395 


± 0.045 


050209 


3.311 


± 


0.274 


0.669 ±0.111 


1.2 ±0.212 


4.24 ± 1.36 


0.663 


± 0.043 


050408 


1.778 


+ 


0.083 


1.103 + 0.005 


0.537 + 0.007 


2.68 + 0.01 


0.261 


+ 0.002 


050509 


1.104 


+ 


0.345 


1.12 + 0.095 


-0.836 + 0.125 


0.64 + 0.28 


0.807 


+ 0.026 


050729 


2.455 


+ 


0.084 


0.172 + 0.059 


1.515 + 0.214 


4.87 + 0.86 


0.447 


+ 0.026 


050807 


2.25: 




0.224 


0.781 ±0.117 


0.15 ±0.117 


1.14 ±0.67 


0.411 


± 0.066 


050922 


2.656 


± 


0.073 


0.285 ± 0.062 


0.742 ±0.132 


4.59 ± 0.08 


0.378 


± 0.007 


051021 


2.386 


+ 


0.067 


1.184 + 0.053 


0.685 + 0.061 


3.04 + 0.65 


0.289 


+ 0.028 


051022 


2.585 


+ 


0.022 


1.995 + 0.012 


1.356 + 0.014 


3.71+0.05 


0.065 


+ 0.001 


051028 


3.103 


+ 


0.152 


0.509 + 0.1 


1.273 + 0.125 


6.72 ± 2.05 


0.522 


+ 0.049 


051211 


2.862 


± 


0.101 


-0.076 ± 0.085 


2.24 ± 0.377 


7.71 ± 1.85 


0.652 


± 0.035 


060115 


2.397 


± 


0.056 


0.91 ±0.025 


1.096 ±0.065 


4.63 ± 0.46 


0.203 


± 0.014 


060121 


2.526 


± 


0.043 


-0.045 ± 0.033 


1.499 ±0.069 


5.37 ± 0.62 


0.253 


±0.018 



